August 4, 2009 



CFHTLS* weak-lensing constraints on the neutrino masses 



I. Terenol'^ C. Schimd^'^, J.-P. Uzan^, M. Kilbinger^, F. H. Vincent^ and L. Fu^'^ 



,2,6 



Argelander-Institut fiir Astronomie, Auf dem Hugel 71, 53121 Bonn, Germany 

Institut d'Astrophysique de Paris, CNRS UMR 7095 & UPMC, 98 bis bd Arago, 75014 Paris, France 
Laboratoire d'Astrophysique de Marseille, 38 rue Joliot-Curie, 13388 Marseille, France 
Universite de Provence - Aix - Marseille I , Marseille, France 

INAF - Osservatorio Astronomico di Capodimonte, via Moiariello 16, 80131 Napoli, Italy 
Shanghai Key Lab for Astrophysics, Shanghai Normal University, 200234 Shanghai, P. R. China. 



Received 3 October 2008 ; accepted 7 February 2009 



Abstract 



Context. Oscillation experiments yield strong evidence that at least some neutrinos are massive. As a hot dark-matter 
component, massive neutrinos shuld modily the expansion history of the Universe as well as the evolution of cosmological 
perturbations, in a different way from cold dark matter or dark energy. 

Aims. We use the latest release of CFHTLS cosmic-shear data to constrain the sum of the masses Yl of neutrinos, 
assuming three degenerate mass states. We also consider a joint analysis including other cosmological observables, 
notably CMB anisotropics, baryonic acoustic oscillations, and distance modulus from Type la supernovae. 
Methods. Combining CAME with a lensing code, we compute the aperture mass variance using a suitable recipe to deal 
with matter perturbations in the non-linear regime. The statistical analysis is performed by sampling an 8-dimensional 
likelihood on a regular grid as well as using the importance sampling technique. 

Results. We obtain the first constraint on neutrino masses based on cosmic-shear data, and combine CFHTLS with 
WMAP, SDSS, 2dFGRS, Gold-set, and SNLS data. The joint analysis yields 0.03 eV < < 0.54 eV at the 95% 

confidence level. The preference for massive neutrinos vanishes when systematics are included. 
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1. Introduction 



The const ruction of a cosmological model (jBondi 19601 : 
lEllis 197 T). must take into account any progress in the un- 
derstanding of the laws of physics. To date, the reference 
model ACDM is based on the standard model of particle 
physics, general relativity, and some additional hypothesis 
about the symmetries of the background spacetime (the 
Copernican principle) to which it is mandatory to add the 
two still unknown components of cold dark matter (CDM) 
and cosmological constant (A). 

We know now that it is imperative to include the 
effects of massive neutrinos, which behave as a warm 
or hot dark matter component depending on their 
mass. Strong evidence of such particles, which have 
not yet been included in the standard model of par- 
ticle physics but who se existance is supported by e.g. 
grand-unified theories (jCalibbi et al. 20 06). emerges from 



experimental results on oscillations of atmospheric, so- 
lar, and accelerator- or reactor-produced neutrinos, such 
as Super-Kamiokande, K2K, M INOS, KamLAND (see 
iGonzalez-Garcia &: Maltoni 20081 . for a review). These re- 
sults indicate a mixing of the three known neutrinos species 
by non-vanishing squared mass differences between the 
mass eigenstatcs and the non- vanishing corresponding mix- 
ing angles. The most recent results are Am^2 '^-^ 8 x 
10~^ eV^, Arn23 ~ 2 x lO^'^ eV^, and mixing angles of 
6*12 ~ 30°, 6*23 45°, and ^13 < 11°. These values yield a 
lower bound on the heavier mass of order 0.06 (0.1) eV for 
normal (inverted) hierarchy, but cannot give the absolute 
mass scale. Forthcoming experiments based on other mech- 
anisms, such as SuperNEMO (conceived to detect the neu- 
trinoless double (3 decay) and KATRIN (designed to probe 
the tritium (3 decay) are designed to measure directly the 
absolute mass of the electron neutrino with a sensitivity of 
0.05 and 0.2 eV, respectively. 
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* Based on observations obtained with 

MegaPrime/MegaCam, a joint project of CFHT and 
CEA/DAPNIA, at the Canada- France-Hawaii Telescope 
(CFHT) which is operated by the National Research Council 
(NRC) of Canada, the Institut National des Sciences de 
rUnivers of the Centre National de la Recherche Scientifique 
(CNRS) of France, and the University of Hawaii. This work is 
based in part on data products produced at TERAPIX and 
the Canadian Astronomy Data Centre as part of the Canada- 
France-Hawaii Telescope Legacy Survey, a collaborative project 
of NRC and CNRS. 



Cosmological observations are compelling in provid- 
ing an independent way to probe the absolute neutrino 
mass. The effect of lig ht massive fermions in cosm ology 
is well understood (see iLesgourgues fc Pastor 20061 for a 
review). The standard cosmological model predicts the ex- 
istence of relic neutrinos that decouple from the primeval 
plasma while ultra-relativistic and produce a cosmic neu- 
trino background, which cannot yet be detected directly. 
Their contribution to the total radiation energy density 
is a fraction 7 /8{Tv/Tj)'^Nes of the contribution of pho- 
tons, and Tj being the temperatures of the neutrino 



and the photon backgrounds, respectively. The effective 
number of neutrinos is constrained to be Ncs — 3.04, if 
there are no extra rela tivistic degrees of fre edom besides 
three active neutrinos ([Mangano et al. 2002f ). and has an 
impact mainly on the primordial nucleosynthesis and the 
time of matt er-radiation equalit y. Analyses of light element 
abundances (ICuoco et al. 2004 } and of the time of equality 
(jKomatsu et al. 2008f ) produce results that agree with this 
value. 

When neutrinos become non-relativistic, their veloci- 
ties begin to decrease from c as T,y(t), due to momentum 
conservation. Assuming that the mass states are degen- 
erate {nil = W2 = ma = m^), the velocity is given by 
Vth = ikBT^(t)/m^c or, since T^(z) — 1.9(1 -I- z) K, by 



va, = 150 (1 + z)(m^/l eV)~^ km/s. 



The transition occurs at 



1 



(1) 



(2) 



For masses smaller than ~ 0.6 eV, the transition occurs 
well into the matter-dominated era and after recombina- 
tion. Therefore, massive neutrinos affect the anisotropics of 
the CMB temperature only through the background evolu- 
tion. However CMB data strongly constrain other cosmo- 
logical parameters that are degenerate with the neutrino 
mass, and is thus invaluable to joint analyses. The result 
would slightly differ for non-degenerate masses, for the nor- 
mal or inverted hierarchi es, but at a negligible le vel for a 
total mass above 0.2 eV (jLesgourgues et al. 20041 ) . 

After the transition, neutrino-density perturbations 
evolve in a similar way to CDM on scales larger than a 
free-streaming length, and are damped on smaller scales. 
The comoving free-streaming length is defined analogously 
to a Jeans length, by replacing the sound of speed by wth, 
i.e., Xfs(t) — 2Tr^/2/Svth{t)/[a(t)H(t)]. It decreases from 
the time of transition onwards as, 

Afs = 2ttx 1.2(1 + zy^^im^/leYy^n^h^y^/^ Mpc, (3) 

which is valid at high redshift. The maximum occurs at the 
transition and is given by, 

/ ID \ — 1/2 

= 27r X 55(an/i^)"^/^ [T^J h-^^pc. (4) 

The density perturbations of the other components, in par- 
ticular CDM, are affected by the presence of massive neu- 
trinos as a consequence of the change in the background 
evolution, the change in the time of equality, and feed- 
back from the neutrino perturbations. Roughly speaking, 
the effect of rrii, = leV neutrinos in the power spec- 
trum is similar to the CDM effect on large scales k = 
27r/A < fcfs(ztr) — 0.01 Mpc~^; a scale-dependent sup- 
pression of the power spectrum amplitude at intermediate 
scales 0.01 Mpc~^ < ^ < ^fs(z = 0) ~ 0.5 Mpc"^ and a 
scale-independent suppression of amplitude at small scales 
k > 0.5Mpc-^ 

Various combinations of cosmological data have been 
used to constrain the neutrino mass. Using CMB data 
alone, the WMAP 5-year analysis obtaied an upper 95% 
confidence limit on the sum of the neutrinos mass es of 
1.3 eV, for a ACDM cosmology (jKomatsu et al. 2008f ). By 
adding galaxy clustering data to CMB data and marginal- 
izing over the galactic bias, the upper limit decreases to 



aroun d 0.7 eV - 1.0 eV, d e pending on the datasets used 
(e.g., iTegmark et al. 20061 : iFogh et al. 20081 ). If instead, 
distance indicators (supernovae and baryon acoustic oscilla- 
tions) were combined with CMB, str onger constraints wer e 
obtained, e.g., < O-^l eV ('Kom atsu et al. 20081 ). 

When all of these probes were used tog ether, the result im- 
proved to around 0.4 eV - 0.5 eV (e.g.. iGoobar et al. 20061 : 
iKristiansen et al. 2007[ ) . Similar results were obtained us- 
ing only CMB and gal axy clustering data bu t assuming that 
the bias was kno wn ([McTavish et al. 2006( 1. or measuring 
it independently (jSeliak et al. 20051 ). Bv applving the lat- 
ter technique to WMAP 5-year data and new galactic bias 
constra ints, an improved value of ^ rrii, < 0.28 eV was ob- 
tained (jDe Bernardis et al. 20081 ). Lyman-a data was also 
added to different combinations of probes, producing re- 
sults of around 0.15 eV - 0.3 eV, in cluding the tighest con- 
straint to date of < 0.17 eV (jSeliak et al. 2006f ). 

In this paper we present the first constraints on neutrino 
masses obtained with cosmic shear data. The possibility 
of using the effect of weak gravitational lensing of back- 
grou nd galaxies by larg e-scale structures, or cosmic shear 
(see iM unshi et al. 20081 for a review), to constrain neu- 
trino properties was first studied in Cooray (1999). Since 
then, various forecasts have been made and future cosmic 
shear surveys in combinations with redshift information 
and CMB data from Planck are expected to measure the 
sum of neutrino masses with an accuracy of 0.05 eV, us- 
ing shear tomographv ( Hannestad et al. 2006[ ) or the full 
3D shear field (Kitching et al. 2008| ). In the present work, 
we use the latest release of the CFHTLS wide survey, 
CFHTLS-T0003, where two-point angular correlations of 
the cosmic shear field were detected on scales ranging from 
1' to 4° (|Fu et al. 20 08). The signal was measured on 
source galaxies with mean redshift Zm = 0.92, implying 
it is produced mainly by the dark matter distribution at 
z ~ 0.4. The measured range therefore probes the non- 
linear matter power spectrum on comoving scales of 0.08 < 
/c < 10 ft. Mpc~^, reaching the quasi-linear regime of largest 
scale corresponding to a physical size of approxiamtely 
1 Gpc. This largest scale probed is still one order of 
magnitude smaller than the largest free-streaming scale, 
/Cfg ~ 0.01 Mpc^^, computed from Eq. ([4]). The data probes 
thus both the scale-independent and the scale-dependent 
suppression effects of the sub-free-streaming regime. 

We proceed by presenting in Sect. 2 the methodology 
followed to compute the spectra. In the same section we in- 
troduce the statistical analysis performed and the datasets 
used. Cosmological constraints are obtained and discussed 
in Sect. 3 and our conclusions are presented in Sect. 4. 



2. Methodology 

2.1. Implementing massive neutrinos in cosmic-shear 2-point 
correlation functions 

We perform a statistical analysis of a mixed dark-matter 
scenario by considering a 8-1-1 dimensional parameter space 
(see Sect. 12. 2p . For each sampled model of this parameter 
space the cosmic-shear power spectrum Pj{£) is computed 
by integrating the non-linear total matter power spectrum 
along the line-of-sight, up to the limiting redshift of 
the survey zum, to account for the distribution of source 




Figure 1. Left: Equation-of-state of fermionic massive particles as function of redshift z for m — 1.0 eV (solid line), 
0.1 eV (dashed) and 0.01 eV (dotted). The transition redshift, Eq. is indicated in each case. The dashed-dotted curve 
shows the lensing efficiency window (arbitrary normalization). Right: Aperture-mass variance (^M^p{6)') with massive and 
massless neutrinos as function of angular scale. The massless neutrinos model (blue, thin lines) is shown in both linear and 
non-linear (halofit) approximations. This indicates that the CFHTLS-T0003 data range (defined by the vertical lines) 
lies mostly in the non-linear regime, where the linear-approximation curve is much lower, and reaches the quasi-linear 
regime. Black, thick curves show three models of massive neutrinos (same masses and line types as in the left panel) with 
identical total mass densities (rim^^ = 0.14). The models for the two lighter neutrino cases are indistinguishable, while 
the case = 1.0 eV produces a significantly lower signal at small scales. 



galaxies n(z). It is given by 



4c4 J, H{z) 
where k ~ £/ Sk{x{^)) £^nd 



il + zrW'{z)P^ik,z), (5) 



W{z) = 



"■"d.'n(.')^-[^(^') 



(6) 



defines the lensing efficiency window function, and S'i<-[x(^)] 
denotes the comoving angular diameter distance as function 
of redshift, which reduces to the comoving radial distance 
x{z) for spatial flatness K = 0. The redshift distribution of 
sources was determined in Fu et al. (2008) from the photo- 
metric redshift estimations of Ilbert et al. (2006). We adopt 
the commonly used fitting function of 

n{z)^N(z/zXexp[-{z/z,)''] (7) 

with best-fit parameters a = 1.47, /3 = 2.15, and Zg = 0.90, 
and a normalization constant N that was determined by 
integrating up to — 6. 

The shear power spectrum is integrated over the ap- 
propriate wi ndow function to obt ain the mass variance in 
apertures e (|Schneider et al. 1998D . 



p^(£). 



(8) 



where J4 is the 4th-order Bessel function of the first kind. 

All quantities relying on the background evolution 
depend on the equation-of-state of the various matter 
components. The neutrino equation-of-state is defined as 
w^{m^,z) = P{'m^,T{z)) / p(m^,T{z)). We compute the 
pressure P and the energy density p at points of a grid 
(toi/, T(2;)), using the fact that the distribution of neutri- 
nos in the phase space is a Fermi-Dirac distribution. The 
result is well approximated by the fitting function 



1 



(1 + z) 0.058 eV 



(9) 



where a = 1.652 and b = 0.561, and is shown in Fig. \T\ (left 
panel). The transition to the non-relativistic regime occurs 
earlier for heavier particles, at the redshift given by Eq. ^ 
and indicated in the figure. The lensing efficiency window 
shows the redshift range probed by our data. At the upper 
end of this range, Wi, is far from zero, especially for low 
masses. There is thus a degeneracy between the character- 
istic redshift of the sources and the neutrino mass, which 
is independent of the well-known redshift-mass degeneracy 
defined by the amplitude of the cosmic shear signal. 

The total matter power spectrum in Eq. ([5]) is computed 
following Hannestad et al. (2006), which is also similar to 
the description ado pted for perturbed quintessence fields 
(ISchimd et al. 20051) . By assuming that the neutrino over- 
densities remain always in the linear regime while CDM and 
baryons grow non-linearly, the total matter power spectrum 
is given by 



Pj"'(/c, z) — 



(10) 

where fi = ili/fl^ with i = {i^, cdm, b} are the density 
fractions of each matter field over the total matter den- 



sity parameter 



rib + ^u- The linear power 



spectra Pj''"(fe,z) are computed using the Boltzmann code 
CAMBQ, to take into account properly the scale-dependent 
growth of the neutrino perturbations. The output is com- 
puted for a grid of redshifts linearly spanning the range 
< 2: < 4 and then spline-interpolated. The non-linear 
spect rum Pr'n'm-i-h is co mputed with the halofit map- 
ping (ISmith et al. 20031) . The Peacock & Dodds mapping 
(jPeacock fc Dodds 1996f ) is unsuitable for this study, be- 
cause it uses a scale-independent growth factor, which can- 
not be defined in the presence of neutrinos. 

The right panel of Fig. [T] shows the aperture mass dis- 
persion in the presence of massive and massless neutrinos, 
computed from the total matter power spectrum, for the 

^ http : 77camb . info] 



same flra- The models are similar on large scales where 
the neutrinos behave like cold dark matter and show a 
scale-dependent deviation on intermediate scales. On small 
scales, which stay always inside the free-streaming length, 
the suppression is constant. Linear perturbation theory 
predicts the small-scale suppression to be, [P^'"(/y,2: = 

0) - Pi!"ff, = 0,z = 0)]/P^^{U,z = 0) 8U 

(|Hu et al. 1998t) . For the non-linear power spectrum, the 
suppression is higher, ar ound —10/^, as predict ed by both 
numerical simulations i Bran dbvge et al. 2008 j and one- 
loop corrections (,Wong 2008, ). For the aperture mass dis- 
persion, the models shown in Fig. [T] have a small-scale sup- 
pression of ~ — 5/jy. 

2.2. Likelihood analysis 
2.2.1. Cosmic shear alone 

The statistical analysis evaluates a mixed dark matter sce- 
nario with nine parameters, including the reduced Hubble 
parameter h, the density of baryons Wb = ^^b^^, cold dark 
matter, Wcdm = ^^cdm^^, and massive neutrinos, uj^, = 
fl^h"^, the effective number of relativistic degrees of free- 
dom A^eff , the optical depth to reionization r, the primordial 
spectral index rig, the rms of the matter perturbations ex- 
trapolated to redshift z — Q and filtered over the 8/i~^ Mpc 
scale, CTgj and the mean redshift of the source galaxies. 
The results will be marginalized over a mean redshift range 
0.78 < Zm < 1.0, corresponding to the 2(t interval of the 
redshift distribution found in Fu et al. (2008). The neu- 
trino physical density parameter relates to the total neu- 
trino mass as Q,^h^ = ^ my/93 eV, and we assume that 
the three neutrino masses are degenerate. 

In a very conservative approach, we shall work in 
a Friedmann-Lemaitre- Robertson- Walker (FLRW) cosmol- 
ogy with cosmological constant as dark energy, with energy 
density fixed assuming spatial flatness. More generally, in 
this framework, dark energy is a fluid-like component that is 
dominant at low redshift and responsible for the recent ac- 
celerated expansion. A quintes sence field, eventually allow- 
ing for non- minimal couphngs (IWetterich 1995HUzan 1999t 
I To echini- Valentin! fc Amendola 2002|) . provides the most 
general alternative that might alleviat e the coinci- 
dence problem (e.g., ICopeland et al. 2006L for a review). 
Alternatively, it is possible to account for the global dy- 
namics without advocating such an additional field, by 
considering back-reac ti on effects of structure formatio n 
(|Ellis fc Buchert 2005t iKolb et al. 200l [Wiltshire 2007D . 
rejecting the Copernican principle ( Uzan et al. 2008[ )~or 
invoking theories of gravitatio n that diffe r from general 
relativity on large scales (see lUzan 20061 for a review). 
Our results do not exclude these possibilities, which re- 
quire ne w, consistent analy ses to be evaluated properly 
(see e.g.. lSchimd et al. 20071 for a weak-lensing analysis of 
quintessence) . 

The log-likelihood is defined as usual to be 

log/: = - J E ^ (Mfp(^.)) {Coy-%A {Mlie,)) , (11) 



where A(Mfp(6'i)) is the difference between the observed 
and the theoretical values of the aperture mass variance 
computed at angles V < Oi < 230' and Cov is the corre- 
sponding covariance matrix. 



We sample the likelihood over a regular grid of the pa- 
rameter space, for a total of ~ 10^ models. The domain 
of the grid is : Wb e [0.0186,0.0249], Wcdm e [0.110,0.152], 
Wy e [0,0.042], h e [0.62,0.80], N^s e [2,4], r e [0,0.225], 
Us G [0.87,1.02], erg G [0.65,0.92]. 

In this work we use the aperture mass vector and co- 
variance matrix obtained in Fu et al. (2008), which was 
measured over 57 square degrees (35 square degrees of ef- 
fective area) in three connected patches of the sky, using 
2.5 X 10^ galaxies with magnitudes 21.5 < zab < 24.5. 
The covariance matrix includes shape noise in the shear 
estimator, non-Gaussian cosmic variance, and the resid- 
ual B-mode. The shear-meas urement pipe line was tested 
on STEP2 simulations (Mass ev" et al. 20071 ) and shown to 
be sl ightly biased, u nderestimating the shear on average by 
2% ()Fu et al. 20081 ). Moreover, this analysis showed that 
highly anisotropic point-spread-functions (PSF) may intro- 
duce a spurious constant shear. The aperture mass statistic, 
being computed with a compensated filter, is unaffect by 
a constant shear, which is the main reason for choosing 
this statistic for the present analysis. In contrast, a con- 
stant shear affects top-hat two-point statistics. The impact 
of redshift-dependent additive and multiplicative shear bi- 
ases, which are present in all KSB-bascd shear estimators, 
on cosmological constraints is studied in Semboloni et al. 
(2008). 

2.2.2. Joint analysis 

In the join t analysis, we a dopt the importance sampling 
technique (jHastings 19701 : iLewis fc Bridle 200"1 ). adding 
CFHTLS-T0003 cosmic shear data to two Monte Carlo 
Markov chains available in the Lambda archive^ for the 
mixed dark-matter scenario. One of the chains was com- 
puted for CMB temperature and polarization anis otropics 
derived from WMAP-5yr (jPunklev et al. 20081 ) data. 
The other chain used a combination of these data 
with b aryon acoustic osci llations from SDSS and 2dF- 
GRS (jPercival et a l. 20071), and Type la supernovae 
from the "Gold" sample (|Riess et al. 2004D and SNLS 
(|Astier et al. 2006D . 

Importance sampling consists of estimating a target dis- 
tribution (a joint posterior distribution in our case) by sam- 
pling an auxiliary distribution. We assume that the public 
Monte Carlo Markov chains provided by the WMAP team 
have enough resolution at the region of parameter space in- 
tersecting the cosmic shear posterior distribution. This is a 
reasonable assumption since they are the dominant contri- 
bution to the joint constraints. They are thus useful auxil- 
iary distributions and are biased distributions of the target 
distributions. To obtain the two unbiased joint distributions 
(WMAP-hCFHTLS and WMAP-fBAO-l-SNe-l-CFHTLS), 
we scale the chains by multiplying the weight of each chain 
element by its likelihood with respect to the cosmic shear 
data. 

This method allows us to accelerate the computation 
compared to sampling directly the joint posterior with a 
grid or a Markov chain, since it requires only the compu- 
tation of the cosmic shear likelihood for each model of the 
chains. Each model of the chains has a constant A'oft = 3.04, 
hence in the joint analysis, we consider only 7+1 physically 
independent parameters, i.e., the remaining 7 cosmologi- 
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cal parameters used in the cosmic shear analysis and the 
lensing sources redshift. 

3. Results 

3.1. CFHTLS-T0003 alone 

We now discuss the analysis of cosmic-shear data alone. 
Figure [2] shows the marginalized confidence contours in 
the {as,^m) plane. The degeneracy direction in the case 
of massless neutrinos is fitted by 

as (f7m/0.25)°-^2 ^ 0.744 ± 0.052. (12) 

This value was found by marginalizing over h but for fixed 
Wb and Us to be more directly comparable to the result of Fu 
et al. (2008) for the same data. Our result is slightly lower, 
but within la. There are several differences between the two 
analyses that account for this difference. In particular, since 
we are interested in combining our results with CMB con- 
straints, we explored a narrower grid in (wb, Wcdm, Wj^), cor- 
responding to a grid interval of Sim G [0.2, 0.6]. The bound- 
ary effects slightly increase the curvature of the contour, 
which is most accurately described with a larger exponent 
power-law. The comparison of both results is particularly 
sensitive to this choice because here the confidence intervals 
are computed as volumes of the posterior distribution con- 
taining a certain fraction of the total probability, whereas 
in Fu et al. (2008) the results are obtained from the val- 
ues of Ax^. The theoretical modelling also differs. Here we 
use Eq. (fTU]) . while in Fu et al. (2008) the power spectrum 
was computed from the Eisenstein & Hu transfer function 
(jEisenstein fc Hu 19981 ). 

In the presence of massive neutrinos, the contours are 
shifted towards the right. This follows from the r^cdm — 
degeneracy, and shows that an increase in fli, is compen- 
sated by an increase in ricdm, confirming that the data are 
in the sub-free-streaming-length regime, where the effect of 
the neutrinos opposes that of the cold dark matter. The 
contours also broaden, since many more preferred models 
exist now. For example, on angular scales of between lO'and 
100', (Mfp(6l)) of the best-fit models of the full grid and the 
massless neutrino sub-grid, differ by less than 5%. As the 
contours broaden, their shape remains unchanged, imply- 
ing that the CDM effect on the cosmic shear signal is more 
significant than that of the neutrinos. The degeneracy di- 
rection in the case of massive neutrinos is now fitted by 

as (flm/0.35)°-^'' = 0.711 ± 0.076 , (13) 

where we chose a higher pivot fim. 

Figure [3] (left panel) shows the massless and massive 
neutrino best-fit models, (m^ ~ eV, fim = 0.30) and 
(m^ = 0.53 eV, flm = 0.44), respectively. They have the 
same values of as and h and differ in rig. The difference 
between the two models, shown in the inset normalized 
by the massless neutrino model, increases on larger scales 
with the approach to the free-streaming-length. This in- 
dicates that combining cosmic shear measurements at sub- 
and super-free-streaming-scales, might break the ricdm ~ 
degeneracy. Shown explicitly in Fig. [3] (middle panel) , this 
degeneracy, as well as the similarly-oriented as — Sljy cor- 
relation (not shown), confirm that an increasing neutrino 
density decreases the cosmic shear signal. However, the tilt- 
ing effect of the neutrino density on the matter power spec- 
trum around the free-streaming scale may mimic the depen- 
dence on rig, especially if the pivot scale k used is similar to 



b 




0.2 0.3 0.4 0.5 0.6 



Figure 2. Confidence contours (68% and 95%) from the 
aperture-mass dispersion between I'and 230', for fl^ — 
(red, smaller) and marginalized over Ui, (blue, larger). Both 
cases are marginalized over the source redshift distribution 
and the remaining cosmological parameters. 

the free-streaming scale ^ O.Ol/i Mpc~^. Therefore, the 
combination of sub- and super-free-streaming data may not 
have enough information to additionally break the Ug — 
degeneracy. This degeneracy is shown in the right panel of 
Fig. El 

To explore the benefit of a broader range of scales, 
we construct a fiducial model consisting of an extension 
of our best-fit massless neutrino model to 6* = 20° (also 
shown in Fig. [31 left panel), where the flat-sky approxi- 
mation remains valid. We computed a new covariance ma- 
trix using the Schneider et al. (2002) approximation with 
the Semboloni et al. (2007) non-Ga ussian correction an d 
WMAP-5yr mean parameter values (|Dunklev et al. 20081 ). 
We assumed the same sky coverage and galaxy density as in 
CFHTLS-T0003. In this way we can evaluate the benefit of 
using large scales independent of a gain due to better statis- 
tics. We obtain constraints similar to the ones obtained for 
the data up to 230', the extension to 9 = 20° only adding a 
small number of independent points. Since the cosmic shear 
signal on large scales becomes very small, better statistics 
are needed in attempting to break the Slcdm — degener- 
acy with cosmic shear data alone. 

Using either CFHTLS-T0003 data or its extension to 
9 ~ 20°, the marginalized upper bound on the neutrino 
mass is, 

^m^, < 3.3eV (95% C.L.). (14) 

This constraint, while very loose when compared with the 
combined constraints mentioned in Sect. 2, is comparable 
to the ones obtained from SDSS or 2dF galaxy redshift 
surveys alone (see e.g., Kristiansen et al. 2007, who find 
X;to^ < 5.2 eV). 

We note that the 95% C.L. contours in Fig. [21 do not 
close inside the range of the parameters probed. This im- 
plies that our grid limits are an effectively strong prior, in 
particular implying fl^ e [0.2,0.6] and as G [0.65,0.92]. In 
other words, the constraint in Eq. (|14p includes a marginal- 
isation over an arbitrary and relatively narrow range of 
parameters values, and is thus optimistic. Furthermore, 
the derived constraint does not take account of degen- 
eracies with dark energy. In more general scenarios with 
non- vanishing curvature or with non-cosmological-constant 
dark energy, both growth of structure and distances, as 
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Figures. Left: Best-fit models found for massive (m^ = 0.53 eV, fim = 0.44; dashed) and massless (m^ = eV, 
rim = 0.30; dotted) neutrino cases. Error bars are from CFHTLS-T0003. The models are extended beyond the data limit 
to illustrate the model separation at large scales. The normalized difference between the models is shown in the inset, 
where the horizontal line marks a 5% model-separation and the vertical line shows the data limit. Middle and Right: 
Confidence contours (68% and 95%) from the aperture-mass variance, in the {uj^, t^^cdm) and (cji/, rig) planes, marginalized 
over the hidden parameters. 




functions of redshift, depend on both the dark-energy den- 
sity parameter Ode and the equation-of-state Wde- To ef- 
fectively break the O^ — Om degeneracy, it will thus be 
necessary to use a measure that is not degenerate in 
Ode and Wdc- One possibility is to exploit the lens ef- 
ficiency at different redshift bins with tomography (e.g., 
lAbazaiian fc Dodelson 20031) . However, even with tomog- 
raphy, some direct neutrino-dark energy degeneracy will 
remain due to the transition epoch from matter to dark- 
energy domination (jKiakotou et al. 2008t ). 

Regarding the effective number of relativistic species, 
cosmic shear is sensitive to it via the change in the matter- 
radiation equality. A higher N^s implies a longer radiation- 
dominated era, and thus a more efficient suppression of 
growth. Hence, the relative amplitude of small and large 
scales changes and the shear power spectrum tilts. The ef- 
fect is however weak and we find an essentially fiat likeli- 
hood in the probed range of iVcff £ [2, 4], meaning that the 
data do not constrain this parameter. 

3.2. Joint analysis 

We now explore the parameter space by adding CFHTLS- 
T0003 cosmic shear data to each of the Monte Carlo Markov 
chains defined in Sect. 2.2.2, aiming to determine the gain 
achieved by including lensing data. 

The main independent piece of information provided by 
cosmic shear is the power spectrum on small scales. As dis- 
cussed above, it has a distinctive scaling with both neutrino 
and CDM densities, while these two parameters for CMB 
are instead weakly correlated, as seen in Fig. 3] (left panel, 
largest contour). This allows us to obtain a narrow, joint 
constraint, approximatively in the same direction found for 
cosmic shear alone in Fig. [3] (middle panel). The direction 
orthogonal to the contour indicates the combination of pa- 
rameters to which cosmic shear is mostly sensitive to. It 
is well approximated by the linear relation Wcdm — and 
the WMAP5 -I- CFHTLS constraint on this combination of 
parameters is, 



The figure-of-merit (FoM) quantifies the error ellipses 
of correlated parameters and was introduced for the 
parameters of the dark energy equation of state 
()Albrecht et al. 20061 ). It is defined as the inverse area of 
the ellipse approximating the p% confidence level contour, 
centered on the contour centroid, i.e., 

FoM-i = Ax'(p)7r ((711^22 - (Tl2f'^, (16) 

where an is the variance of the parameter i. We consider 
2a contours, for which Ax^(95.4%) = 6.17, and obtain 
FoM ~ 19 for the WMAP5 contour and FoM ~ 44 for 
the combined (WMAP5 -I- CFHTLS) case, corresponding 
to a gain of a factor of 2.3. The two parameters ujcdm and lOh 
arc, however, strongly correlated, exhibiting a narrow, long 
combined contour and the effective gain in the lo^ variance 
is of only 1.2. 

The introduction of distance measurements 
(BAO+SNe) is useful, since its combination with WMAP5 
data (also shown in Fig. HI left panel) defines a narrow 
region in the iuicdun^v) plane that complements that 
of the WMAP5-hCFHTLS data. Both cases thus probe 
orthogonal combinations of the two parameters and, 
furthermore, in a way that is independent of the spectral 
index, since this is already well determined by the CMB 
data alone (as shown in Fig. [H middle panel). Their 
combination thus has the same properties we were looking 
for when discussing the possibility of using sub- and super- 
free-streaming scales from cosmic shear alone, along with 
the need for an independent measure of rig. This means 
that the combination of the two cases (WMAP5 -f BAO -I- 
SNe -I- CFHTLS) breaks the remaining degeneracy, as is 
clear from Fig. |4] (left panel, smaller contour) and from the 
fact that the combined contour has very small correlation 
(ai2 < 0.1%, for un-normalised parameters ujc and u>^). 

Figure [4] (right panel) shows the one-dimensional 
marginalized probability distributions of the neutrino 
masses in the four cases. The upper 95% confidence lev- 
els are Y^m^ < 1.32 eV (WMAP5), J^m^. < 1.22 cV 
(WMAP5 + CFHTLS) and J^m^ < 0-62 eV (WMAP5 + 
BAO -I- SNe). For the fuU combination (WMAP5 + BAO 
+ SNe + CFHTLS), we find that 



Wcdm = 0.106 ±0.006 (Icr). (15) 



0.03eV < ^TO^ < 0.54eV (95% C.L.), (17) 




Figure 4. Joint analysis. Left and Middle: Contours at 95% C.L. of ((jjj,,Wcdm) and (wi^jUs) for the four cases consid- 
ered : WMAP5 afone (largest contour, red), WMAP5-I-CFHTLS (second largest, yellow), WMAP5+BA0+SNe (second 
smallest, green), WMAPS+BAO+SNe+CFHTLS (smallest contour, blue). Right: One-dimensional marginals for J^''^'^ 
showing the 95% C.L with vertical hues, for the same four cases : WMAP5 (dash-dot, red); WMAPS-f^CFHTLS (dash, 
pink) ; WMAP5-hBA0+SNe (long dash, green); and WMAP5+BA0+SNe+CFHTLS (solid, blue). For the latter a lower 
bound was also found). 



with mean 0.31 eV. Interestingly, we obtain a lower bound, 
and thus a preference for massive neutrinos at the 2a level. 

It is worth to point out that the lower bound on ^ 
is strongly dependent on the position of the WMAP5 
-|- CFHTLS contour in the {uj^,ujcdm) plane (see Fig. HI 
left panel). Any systematic effect that underestimates the 
cosmic-shear signal would shift the {finncrs) contour to- 
wards the bottom-left corner of Fig. [5] with a consequent 
shift in the (iOi^^ujcdm) contour towards the bottom-right 
corner of Fig. [1 (left panel). The WMAP5 -I- BAG + SNe 
+ CFHTLS contour would thus descend along the WMAP5 
-|- BAO -I- SNe one producing a positive lower bound for 
the neutrino mass. A similar effect appeared in Allen et 
al. (2003), which also showed a preference for a non-zero 
neutrino mass from a combination of CMB, galaxy cluster- 
ing and X-ray cluster data; that result was in fact caused 
by an underestimation of as, as pointed out by Seljak et 
al. (2005). Similarly, the joint constraint depends crucially 
on BAO data. Komatsu et al. (2008) pointed out a tension 
between SNe and BAO results when the joint SDSS -I- 2dF- 
GRS BAO sample is used. If BAO data overestimates the 
total matter density, then an unbiased WMAP5 -I- BAO -I- 
SNe contour in the {uj^, ojcdm) plane (Fig. HI left panel) will 
be shifted towards the bottom-left corner, weakening the 
preference for massive neutrinos. 

3.3. Robustness of the constraints 

We check here the robustness of our results against system- 
atics of the cosmic shear data, that have not been included 
in the analysis so far. For this purpose, we first assume the 
data vector is biased as a consequence of an underestima- 
tion in the shear, as indicated by the STEP simulations. 
We assume a 2% redshift-independent bias in the shape 
measurements, which translates into a scale-independent 
4% error in the data vector, and repeat the joint analy- 
sis. The resulting combined contour shifts upwards along 
the WMAP5-t-BA0-|-SNe degeneracy in agreement with the 
discussion in Sect. 3.2. As seen in Fig. [S] (left panel, sohd 
lines) , this small bias is enough for the evidence of massless 
neutrinos to be lost and we obtain at the 95% C.L. : 

^m^<0.53eT/. (18) 



Other sources of errors that can affect cosmic-shear 
cosmological constraints are: contamination from intrin- 
sic alignments, in p articular shear-shape correlations 
()Hirata fc Seliak 2004 ) for no precise theoretical modeling 
yet exists; uncertainties in the photometric redshifts; and 
extra PSF residuals not predicted by the STEP simulations. 
Modelling of some of the CFHTLS-T0003 systematics is in- 
cluded in the analysis of Kilbingcr et al. (2008). Here we 
consider the rougher approach of marginalizing over a scale- 
independent calibration factor, accounting for all possible 
sources of biases, with the goal of finding a threshold of con- 
tamination above which the CFHTLS data do not improve 
the combined constraints. 

We find a threshold of 25% for the aperture mass dis- 
persion data. The corresponding combined contour in the 
(uj^, ujcdm) plane is shown in Fig. [5] (left panel, dashed 
lines). The contour is now elongated and similar to the 
WMAP-l-BAO-hSNe one of Fig. 11 (left panel). If the data 
contains systematics of amplitude smaller than 25% of the 
signal, we can quote conservative combined constraints by 
marginalizing over the amplitude of the systematics. The 
final 95% C.L. for the joint analysis will be intermedi- 
ate between the most optimistic case (no bias) of 0.03 eV 
< J2''^f < 0-54 eV and the worst case (no contribution 
from lensing) of ^rrii, < 0.62 eV. For example, in the 
case of systematics of 10%, the conservative constraints are 
^rrXj/ < 0.58 eV. These marginalized constraints, not in- 
cluding systematics with a scale- or a redshift-dependence, 
are approximative since massive neutrinos affect both the 
amplitude and the shape of the cosmic shear correlations. 
Only on smaller scales (< 10' , see e.g. Fig. [1] right panel) 
is the suppression scale-independent. 

We can also assess the robustness of the result by 
comparing it with a forecast. For this, we compute a 
systematics-free cosmic-shear covariance matrix for our sur- 
vey size using the results of Schneider et al. (2002) and 
Semboloni et al. (2007). This covariance matrix differs from 
the one we have used so far, which was computed by tak- 
ing into account the true galaxy positions and weights. The 
results do not change when we replace the covariance matri- 
ces, as shown in Fig. [S] (middle and right panels, solid lines). 
To ensure that the forecast is completely independent of the 




Figures. Tests of the robustness of the resuh. Marginahzed joint constraints (68% and 95% C.L.) obtained using 
the actual data (filled contours) and various assumptions (open contours). Left: Impact of systematics, assuming a 4% 
underestimation (solid) or marginalizing over a ±25% calibration bias (dashed). Middle and Right: Forecasts using a 
synthetic covariance matrix and the CFHTLS vector data (solid) or the same covariance matrix with a WMAP5 fiducial 
model (dashed). 



data, we redo it by further replacing the data vector by a 
fiducial model (the WMAP5 mean, iDunklev et al. 2008^ . 
This time wc obtain a shift in the contour location (Fig. [51 
middle and right panels, dashed lines). Again the evidence 
for massless neutrinos is lost, this time due to the higher 
Wcdm value of the WMAP5 fiducial model. This behaviour 
mimicks a correction for an eventual overestimation of the 
source redshifts. Indeed, if the effective redshift is lower, 
the models will have a higher value of fim for the same 
amplitude as, shifting the combined Wcdm — oJi, contour up. 
It likewise mimicks a correction for an underestimation of 
the shear measurements. The combined constraint on the 
neutrino mass, when using the WMAP5 fiducial model, 
strengthens to '^m,, < 0.44 eV. We note this is also the 
result that would be obtained if the data were corrected for 
a negative calibration bias of 20%. 

In summary, the possible systematics do not change the 
angle between the cjcdm — ^^i^ degeneracy in the WMAP5 + 
CFHTLS and WMAP5 + BAG + SNe cases, only shift or 
broaden the combined contour position and size. Eventual 
corrections for an underestimation of the data vector or 
for an overestimation of the source redshifts would produce 
both a tighter upper bound and a looser lower bound on 
the neutrino mass. 

The matter power spectrum also contributes as a source 
of uncertainties. The prescription used, Eq. (fTU]) . does 
not consider clustering of neutrinos on CDM structures, 
which occurs when the neutrino thermal velocity drops be- 
low the velocity dispersion, v, of forming clusters; for in- 
stance, it takes place at z ~ 2.3 for v ~ 1000 km/s and 
rriv ~ 0.5 eV. The corresponding neutrino halo profile is 
flatter in the ce ntre than a pu re CDM Navarro-Franck- 
Whitc ( Ringwald fc Wong 200l ). Including it in the 1-halo 
term of the halo model, Abazajian et al. (2005) showed 
that it decreases the non-linear matter power spectrum. 
Accordingly, in this work the term ^'"dm+b expected to be 
overestimated by ~ 1% on scales around k = 0.5 ft. Mpc~^ 
for ~ 0.5 eV. Alter natively, res ults from perturbation 
theory with neutrinos (jWong 2008t ) also indicate an over- 
estimation of the matter power spectrum, as already men- 
tioned in Sect. 2.1. 

The non-linear power spectrum is computed with the 
HALOFIT, which declares a 3% uncertainty on scales k < 



10 h/Mpc at z < 3 (jSmith et al. 2003D . which are thus 
of the same order as the supposed shear measurement 
bias. This uncertainty is effectively both redshift- and 
scale- independent over the data redshift and scale ranges. 
Marginalizing over a 5% uncertainty, the corresponding 
joint constraint is '^mi, < 0.56 eV. Furthermore, the 
HALOFIT does not take account of the effects of cooling 
baryons and hot intra-cluster baryons, which are degenerate 
with the neutrinos, affect ing the pow e r spectrum by a few 
percent on sma ll scales (jWhite 20051: IZhan fc Knox 2004t 
iJing et al. 20061 ). 

4. Conclusions 

We have investigated the potential of cosmic shear to con- 
strain the mass of neutrinos. In the sub-free-streaming 
regime, the constraining power originates, for a fixed den- 
sity of baryons, in the tendency of relativistic (hot) dark 
matter particles to escape from collapsed regions. Therefore 
additional amounts of CDM are needed to produce the 
same cosmic-shear distortion, shifting the f2ni — erg de- 
generacy towards larger values of i^m with respect to the 
analysis without massive neutrinos (Fig. [2]), producing a 
(^cdm — degeneracy favoring higher amounts of CDM for 
higher amounts of massive neutrinos (Fig. [3]) . The analysis 
of CFHTLS-T0003 data alone yields a loose constraint of 
^m,y < 3.3 eV at the 95% C.L., for our particular choice 
of priors. 

We have explored larger angular scales using a synthetic 
data vector extended to 20° to explore the possibility of 
breaking the Wcdm — degeneracy using only cosmic shear. 
This analysis did not predict an improvement in the results, 
which would require higher signal-to-noise ratio at large 
scales. 

The introduction of CFHTLS-T0003 data in a 
WMAP5-|-BA0-f SNe analysis provides an interesting com- 
bination, breaking the Wcdm — degeneracy present 
in that analysis (Fig. [4] , left panel). The joint anal- 
ysis yields 0.03 eV < Y.^^ < 0.54 eV at the 95% 
confidence level, marginally excluding massless neutrinos 
(Fig. m right panel). The preference for massive neutri- 
nos is lost when the data is corrected for a possible un- 
derestimation of the shear signal, which simulations in- 



dicate may be around 4%. In this case, the result shifts 
to < 0.53 eV. The CFHTLS data also shows 

some disagreement with the result of WMAP5. Indeed, 
if a cosmic-shear fiducial WMAP5 model was used, the 
expected combined WMAP5+BA0+SNe+CFHTLS con- 
straint would be tighter: ^ nii, < 0.44 eV. Finally, we 
showed that CFHTLS-T0003 data do not improve the 
WMAP5-|-BA0-|-SNe constraints on the neutrino mass if 
they contain a bias larger than 25%. 

After submission of the manuscript, a similar study 
was submitted by Ichiki et al. (2008), also using CFHTLS- 
T0003 cosmic shear data to constrain the neutrino mass. 
Both analyses take a similar approach, the main differ- 
ences being: they use the cosmic shear correlation func- 
tion, which is more affected by systematics, as opposed 
to the aperture-mass dispersion; the non-Gaussian shear 
covariance matrix is also computed in different ways, as 
is the likelihood analysis. Their neutrino mass constraints 
from cosmic shear only are weaker, which is consistent with 
the fact they impose weaker priors, and the results for 
a WMAP5-|-BA0-f SNe-t-CFHTLS joint analysis are very 
similar. 

Besides the neutrino mass, we also considered the effec- 
tive number of relativistic degrees of freedom, Ncs- CMB 
anisotropics are by far more sensitive than cosmic shear 
and CFHTLS-T0003 data produced an almost flat likeli- 
hood with respect to this parameter. This prevented us 
from using the as — A^eff degeneracy, caused by a higher 
A^off delaying the matter-radiation equivalence and the cor- 
responding defreezing of matter perturbations, which have 
less time to grow. 

In summary, the CFHTLS cosmic-shear has already suf- 
ficient statistical precision for the accuracy of the results to 
be affected by systematics. The CFHTLS lensing system- 
atics collaboration is currently undertaking a very detailed 
and lengthy analysis of those systematics. The statistical 
precision is however still insufficient to allow us to improve 
current constraints on the neutrino mass. Future ground- 
based and space-borne observations (such as KIDS, Pan- 
STARR S PES, LSST, SNAP -L, JDEM, or EUCLID sur- 
veys; see lPeacock et al. 20061 ). with better statistics, larger 
scales, and also the use of tomographic techniques, will al- 
low to perform more elaborate analyses, for instance en- 
abling to relax the assumption of degeneracy between mass 
eigenstates. 

From the theoretical point-of-view, a more suitable com- 
putation of the matter power spectrum in the non-linear 
regime will be mandatory, to take more careful account of 
the effects of massive neutrinos without relying on map- 
pings based on ACDM A^-body simulations. This might 
be a chievable using e ither standard perturbation theory 
Ce.g..[Saito et al. 20081) or a resummation s cheme (see e.g., 
ICrocce fc Scoccimarro 20061 : iPietroni 20081 ). It will be in- 
teresting, eventually, to investigate the degeneracy with the 
low-redshift physics (z < 20) , such as reionization and dark 
energy, relevant if < 0.01 eV. On this mass scale, the 
normal/inverted hierarchy of mass states can also be dis- 
tinguished. More elaborate models allowing for an inter- 
action between massiv e neutrinos and a quintessence field 
(|Brookfield et al. 20061 ). or in which mass- varying neutrinos 
behave as a negative pressure flui d, which could then be the 
origin of the cosmic acceleration ()Fardon et al. 200^ , may 
then also be considered. 
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